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The ability to engineer both linear and non-linear coupling with a mechanical resonator is an 
important goal for the preparation and investigation of macroscopic mechanical quantum behav- 
ior. In this work, a measurement based scheme is presented where linear or square mechanical 
displacement coupling can be achieved using the optomechanical interaction linearly proportional 
to the mechanical position. The resulting square displacement measurement strength is compared 
to that attainable in the dispersive case using the direct interaction to the mechanical displace- 
ment squared. An experimental protocol and parameter set are discussed for the generation and 
observation of non-Gaussian states of motion of the mechanical element. 



I. INTRODUCTION 

Currently, the main approaches to cavity optomechan- 
ics p] can be divided into two categories - 'reflective' and 
'dispersive'. In each approach the mechanical and opti- 
cal degrees of freedom are coupled via radiation pressure 
and the dependence of the cavity resonance frequency on 
the mechanical position. The first approach is depicted 
in Fig. [lja), where the optical field is reflected from a 
mechanical element and the change in cavity frequency 
and hence interaction Hamiltonian are linearly propor- 
tional to the mechanical position. Optomechanical real- 
izations of this approach include deformable Fabry-Perot 
cavities and optical whispering gallery mode resonators, 
which are discussed in Ref. pQ. The second approach 
is depicted in Fig. [ljb), here a mechanical element is 
positioned within an optical field and partial reflection 
from both sides gives rise to a dispersive interaction. 
In this arrangement the cavity frequency varies periodi- 
cally with mechanical displacement. This can be used to 
give a linear or quadratic position dependent cavity fre- 
quency change if the mechanical element is positioned at 
an anti-node or node of the field respectively. The ability 
to select between linear or quadratic displacement cou- 
pling provides considerable versatility and thus 'disper- 
sive' optomechanics is an exciting candidate to observe 
and explore quantum mechanical phenomena of macro- 
scopic resonators. Optomechanical realizations of this 
approach utilize a dielectric membrane [5] or trapped cold 
atoms ;3] , positioned within an optical cavity and experi- 
mental work is underway to realize this with an optically 
trapped microsphere 0]. The quadratic mechanical po- 
sition coupling offered by dispersive optomechanics pro- 
vides a route to observe quantization in mechanical en- 
ergy [2]. Moreover, such quadratic coupling can also be 
used for cooling and squeezing of the mechanical element 
[5] and it can be enhanced by using additional optical 
spatial modes, which also even allows quartic interaction 

In this paper, a scheme is presented which uses the 
optomechanical interaction linearly proportional to the 
mechanical position where, despite the form of the in- 
teraction, effective couplings to displacement and dis- 



placement squared are achieved. Here optical pulses 
that are short compared to a mechanical period are used 
and the square displacement coupling is obtained by ex- 
ploiting the non-linear optical dependence of the inter- 
action. This interaction has been linearized in much of 
the present literature, but continuous non-linear optome- 
chanical interaction has recently been studied resulting 
in non-classical states of light [7] and of the mechani- 
cal oscillator [3]. Also, in the spin ensemble commu- 
nity, working beyond the linear regime has been used 
for non-Gaussian quantum state preparation 9 . The 
optomechanical setup considered here is shown in Fig. 
[ljc), where an optical pulse in a coherent state interacts 
with an optomechanical system and is then measured via 
homodyne detection. Following the interaction, Wigner 
reconstruction of the optical subsystem of the optome- 
chanical entangled state, would yield a 'scimitar state' 
shown in Fig. [ljd). The form of this state can be read- 
ily understood as the mechanical position fluctuations 
(including quantum fluctuations) rotate the optical field. 
For small rotations, one sees from Fig. Qd) that mea- 
surement of the optical phase quadrature allows for a 
measurement of mechanical position, however, of partic- 
ular interest here, measurement of the amplitude quadra- 
ture may give outcomes which could have resulted from 
two distinct mechanical positions. This is due to an ef- 
fective displacement squared coupling, which can be used 
for non-Gaussian state preparation. In Ref. jlOj it was 
discussed how measurement of the optical phase quadra- 
ture can be used to perform quantum state tomography 
of the motional state of the mechanical resonator and 
generate conditional mechanical squeezed states. Thus, 
being able to select between displacement and displace- 
ment squared measurements provides the tools to gen- 
erate non-Gaussian mechanical quantum states and per- 
form state reconstruction simply by choosing the phase 
in the homodyne interferometer as is shown in Fig. |lje) . 



II. MODEL 

The optomechanical Hamiltonian with linear mechan- 
ical position coupling in the optical rotating frame at the 
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FIG. 1: Cavity optomechanics is currently realized using 
'reflective' (a) and 'dispersive' (b) approaches. The interac- 
tion in the former is proportional to mechanical displacement 
Xm, however, in the latter the interaction can be tuned to 
being proportional to Xm or X M . Pulses of light may be 
used to probe and manipulate the motional state of the me- 
chanical resonator. Considered here is a pulse incident upon 
an optomechanical system with an interaction proportional to 
Xm and then an optical quadrature measurement performed 
via homodyne detection (c). Following the interaction, the 
Wigner function of the optical field (d) is scimitar shaped due 
to mechanical position induced optical rotations. The distri- 
bution of the initial coherent state is indicated by the dashed 
circle. (The parameters for this plot were chosen to exagger- 
ate the curvature.) For small rotations, the phase quadra- 
ture Qp is proportional to Xm and the amplitude quadrature 
carries X M information. The amplitude quadrature measure- 
ment outcome indicated by the red dashed line may have orig- 
inated from two distinct mechanical positions, which provides 
a means for superposition preparation. By choosing the phase 
in homodyne detection (e) one can use the amplitude quadra- 
ture for X M measurements and quantum state preparation 
(QSP) or use the phase quadrature for Xm measurements to 
perform quantum state tomography (QST). 



= u> M ^ b - g Vln V2a^ aX M -iy / 2KN p a in (a-a^), (1) 



cavity frequency including coherent resonant drive is 

Him 

where the optomechanical coupling rate, which is realiza- 
tion dependent, is of the form gy ln — ujlXq/L. The cav- 
ity field's resonance frequency, annihilation operator and 
amplitude decay rate are ujl , & and n respectively, and 
L is the cavity length. The mechanical zero-point exten- 
sion is xq = y/h/2mu)M, where tJM,b,m and Xm(Pm) 



are the mechanical element's eigenfrequency, annihila- 
tion operator, effective mass and position (momentum) 
quadrature, respectively, where a single mechanical mode 
is considered. The input pulse has mean photon number 
N p and is described by a; n , the normalized envelope, i.e. 
Jdtaf n (t) = 1, which is assumed real. 

During the interaction, which is short with respect to a 
mechanical period, i.e. k 3> the mechanical position 
is considered constant and the optical and mechanical 
equations of motion can be solved independently of one 
another. Immediately after the pulse interaction the me- 
chanical position is unchanged, i.e. Xffi = X M , however 
optomechanical entanglement is generated and correla- 
tions are established between the mechanical momentum 
and the optical intensity, Pffl = Pfi + v2<?iin Jdta^a. 

The intra-cavity field evolves during the non-linear op- 
tomechanical interaction according to 
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where the field is rotated in proportion to the mechanical 
position. This can be immediately solved exactly [llj . 
however, in this work the solution is approximated as 
the rotation is assumed small and the mean of the field 
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where the dimensionless temporal mode functions 0:0,1,2 
are introduced [12] . The phase quadrature of the intra- 
cavity field contains information on the mechanical dis- 
placement and the amplitude quadrature carries informa- 
tion of the mechanical displacement squared, Fig. |lje) . 
Measurement of these quadratures can be performed 
by time-domain homodyne detection of the output field 
a out = — a in . Homodyning the amplitude quadra- 

ture is described by Qx = \/2 Jdt aLo(t)X1 ut (t), where 
X£ ut = 2 -1 / 2 (a out + aj ut ) (similarly Qp describes phase 
quadrature detection). For an optimal measurement of 
X M , (Xm) one chooses the local oscillator pulse a^o 
to have an amplitude directly proportional to CV2, (c*i)- 
The mean of the amplitude quadrature measurement is 
(Qx) = Qx ~ Xx(X M ), where the first term is the 
contribution from cvo and Xx is the square displace- 
ment measurement strength. For convenience, the ho- 
modyne measurement outcome is re-written as AQx = 
Qx ~ Qx- The optimal single pulsed measurement of 
Xm is achieved with an input drive with a Lorentzian 
spectrum which matches the natural decay of the cavity 
|10) . The square displacement measurement strength is 
optimal when af n (uj) = (37t) _1 8k 5 /(k 2 +u> 2 ) 3 , which is 
not Lorentzian due to the higher order nature of the inter- 
action considered here. This gives Xx — y / 42A^ p gf in / ' n 2 . 

This kind of pulsed interaction and measurement is 
well suited to being described with the use of measure- 
ment operators as outcome probabilities and conditional 
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mechanical states can be readily determined [T3j . Homo- 
dyne detection of the amplitude quadrature has the out- 
come probability density Pr(AQx) = T^mC^x^xPm)^ 
where is the corresponding measurement operator. 
In this pulsed regime AQx has mechanical dependence 
only on Xm which allows T X T X to be interpreted as an 
outcome probability density conditioned on a mechanical 
position. For the coherent optical drive considered here 
one obtains 



T X (X M ,AQ X ) 



n -i/4 e *n lin x M exp 



-(AQx - XxXm? 



(4) 



where the mean momentum transfer is Oii n 
{5V2/3)N p g nn / K . 



III. COMPARISON TO THE DISPERSIVE 
QUADRATIC INTERACTION 

Before proceeding to a discussion of the mechani- 
cal states of motion that can be prepared with T x , 
the square displacement measurement scheme introduced 
above is compared with the dispersive case. The Hamil- 
tonian from Ref. [5] for optomechanical systems with a 
dispersive element positioned so that the cavity frequency 
varies quadratically with the position of the element, in 
the optical rotating frame at resonance, including drive 
is 



ui M tfb + g^aXlj - i^2nN p a in (a - a^), (5) 



where the quadratic optomechanical coupling rate is 
,g sq = (16w 2 cxl/LX 2 ) v / 2(l - r), r is the (field) reflectiv- 
ity of the dispersive element and A is the optical wave- 
length. The phase quadrature of an optical pulse inci- 
dent upon such an optomechanical system will be dis- 
placed in proportion to X M and it is readily shown that 
for a homodyne measurement of the phase quadrature 
with outcome Qp the measurement operator is T sq = 
w -i/A e -iQ B<l xSX M ex p[-i (Q P + XsqXi;) 2 }, which has 
recently been used in Ref. [13]. After pulse shape op- 



timization, Q sq — 3N p g sq /n and Xs 



^10N p g sq /n. 



Comparing the measurement strengths for the dispersive 
direct X M interaction and the effective X M from the lin- 
ear interaction for identical N p and A gives 



Xx 
Xsq 



7T T f 



x 2 q y/2(l-r) 



(6) 



where the cavity finesses and mechanical zero-point ex- 
tensions are distinguished by subscripts for the two op- 
tomechanical cases. Remarkably, using the optomechan- 
ical interaction which is linearly proportional to Xm and 
optical amplitude quadrature measurements gives an X M 
measurement strength which scales more favorably than 



that available with the direct X M interaction in disper- 
sive optomechanics. This, in combination with the mea- 
surement based selectability between linear or quadratic 
couplings offered here are the main results of this work. 



IV. EXPERIMENTAL PROTOCOL AND 
DISCUSSION 

Jacobs and colleagues discussed the preparation of su- 
perposition of the position of a mechanical resonator via 
X M measurements [15], work which has also recently 
been extended to include feedback control of the superpo- 
sition separation |16) . Such benchmark quantum states 
show striking differences between classical and quantum 
behavior and are thus highly sort experimentally to study 
the quantum mechanical properties of macroscopic ob- 
jects [17HI9] . In the following, an experimental protocol 
and a parameter set are discussed to prepare and observe 
spatial superposition of a massive mechanical resonator 
using the non- linear interaction and measurement T^. 
A measurement on a variety of experimentally accessible 
initial states is considered and the resulting conditional 
and unconditional mechanical states of motion arc deter- 
mined. 

As the spectrum of measurement outcomes is continu- 
ous it is not experimentally possible to post-select from 
many experimental runs on a single outcome, rather a 
window must be used. The mechanical state conditioned 
on outcomes occurring in the window AQx ± w/2 (la- 
beled by w) is 

p * = pTR h AQ ' x T x(AQJc)pS4(AQx)> (?) 

where Pt(w) — fdAQx Pr(AQx) is the probability of 
obtaining an outcome in the window. The mean mea- 
surement outcome for a mechanical thermal state with 
thermal occupation n is (AQx) = Xx (1/2 + n). As the 
mean is greater than zero some insight is gained into the 
form of Pr(AQx), which is a non-Gaussian function with 
a large wing for positive outcomes which increases for 
larger mechanical position variance. 

In Fig. the action of T x is considered on three me- 
chanical Gaussian states: the ground-state, a thermal 
state and a momentum squeezed state. One may sus- 
pect that quite a narrow window for conditioning must 
be used in order to achieve significant coherence between 
the superposition components, however, conditional me- 
chanical states, prepared from high-purity initial states, 
show strong quantum coherence even for relatively large 
conditioning windows. For example, the conditional me- 
chanical state shown in Fig. [2^b) exhibits strong Wigner 
negativity even for w = 0.8, which allows the use of 15% 
of the measurement outcomes. This plot also reveals the 
interesting feature that the negative regions are 'curled 
around' positive regions, a feature which is not seen in the 
more commonly studied superposition of coherent states. 
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TABLE I: An experimentally accessible set of parameters to 
achieve unity square displacement measurement strength. 



Optical wavelength: 


A 


1064 


[nm] 


Mechanical effective mass: 


m 


40 


[ng] 


Mechanical eigenfrequency: 


ujm/2ty 


2 


[kHz] 


Cavity finesse: 


T 


5 x 10 4 




Photon number per pulse: 


N p 


1.7 x 10 9 




Cavity length: 


L 


750 


[Mm] 


Mechanical ground-state size: 


x 


10 


H 


Optomechanical coupling: 


.giin/27T 


3.8 


[kHz] 


Single photon strength: 




1.9 x 10" 3 




Quadratic pos. meas. strength: 




1.0 




Separation (n = 0, AQx = 1.5): 


S 


2.0 





This arises due to the population components having an 
asymmetric distribution about their peaks, specifically, 
there is a broader wing nearer to Xm — and a sharper 
edge on the other side. This form of the population com- 
ponents is more clearly seen in Fig. [2^e), which is the 
conditional mechanical state starting from a low occu- 
pation thermal state. When the population components 
have a more symmetric Xm distribution about their peak 
the interferences no longer curl as strongly, as is seen in 
Fig. [2](li) the conditional state starting from a squeezed 
state. 

Since a measurement of the optical amplitude quadra- 
ture erases all the linear displacement information gained 
during the interaction and (AQx) > 0, the unconditional 
(i.e. all measurement outcomes arc ignored) mechanical 
state = dAQx T x p M T x , is also non-Gaussian, 
however, mixed. 

The superposition separation 6 is defined as the dis- 
tance between the maxima of the two population com- 
ponents. This depends on the initial mechanical distri- 
bution, the measurement outcome and the square dis- 
placement measurement strength. For Gaussian initial 
mechanical states with a standard deviation a in their 
position spread, the superposition separation is 

Experimental progress in optomechanics is steadily ap- 
proaching the regime where the important parameter 
5iin/K, which quantifies the mechanical momentum dis- 
placement by a single photon (for k 3> ^>m) approaches 
unity. In this work Xx scales with the square of this pa- 
rameter and the pulsed position measurement strength 
for mechanical quantum state tomography |10j scales lin- 
early with this parameter. In present-day experiments 
however [20] . gn n /K <C 1, which this work overcomes by 
utilizing large coherent amplitudes in order to achieve 
sufficient coupling to prepare and observe non-Gaussian 
mechanical states of motion. To ensure short interaction 
the cavity decay rate is chosen as k = 10 3 ujm, which 




FIG. 2: Mechanical Wigner functions of initial states (left), 
conditional states (center) and unconditional states (right). 
(The plot range is ±5 for all axes. Color scale: black repre- 
sents zero magnitude, blue for positive values and red for neg- 
ative values.) The initial states are the ground-state, n = 
(a), a thermal state with n — 2 (d) and a momentum squeezed 
vacuum state with squeezing parameter r = 0.5 (g). Condi- 
tional states prepared with Tx acting on the corresponding 
initial states with \x = 1, AQx = 1.5, w = 0.8 are shown in 
(b,e) and AQx = 6.4 has been used in (h). The probabilities 
of obtaining an outcome in the windows used above are: (b) 
14.9 %, (e) 14.5 %, (h) 1.1 %. Note the disappearance of neg- 
ativity - a quantum to classical transition - for initial thermal 
occupation (e) and if the measurement outcomes are ignored 
(c,f,h). 



for a desired finesse sets the cavity length required. In 
Table [I] a list of parameters is provided for a deformable 
Fabry-Perot optomechanical system with a kHz scale me- 
chanical resonator. 

The protocol for quantum state preparation and quan- 
tum state tomography comprises three steps: i) an ini- 
tialization stage of mechanical pre-cooling and/or squeez- 
ing. Since k ^> ujm is required here and low frequency 
mechanical resonators are considered active-feedback 
cooling is most suitable [3TJ [S3] . Alternatively, in this 
regime squeezing and purification can be achieved with 
the use of conditional measurement |10j . Additionally, 
squeezing can be achieved by applying a parametric mod- 
ulation to the mechanical device 23 . ii) Following this, 
a pulse is injected into the optomechanical cavity to 
realize T x and the measurement outcome is recorded. 
At this point, the mechanical oscillator has gained the 
momentum f2r m which after one quarter of a period of 
free evolution shifts the cavity resonance frequency by 
Aojq = giinV^^im- As this can be much larger than k 
any subsequent pulse will not resonantly drive the cav- 
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ity. In order to overcome this, a two pulse preparation 
sequence can be used where a second pulse follows after 
half a mechanical period of free evolution to cancel the 
mean momentum gained by the resonator. In this case 
one applies twice where both outcomes are recorded, 
thus strengthening the measurement of X^ 4 . This pro- 
cedure requires a good degree of optical amplitude sta- 
bility, which is anyway necessary for T x measurements. 
During the free evolution the appropriate master equa- 
tion is solved to determine the mechanical state immedi- 
ately prior to the second measurement. However, as dis- 
cussed below, given the parameters considered here the 
mechanical bath coupling is not expected to play a strong 
role during this timescale. in) With the resonator state 
near the origin of phase-space, quantum state tomogra- 
phy, as discussed in Ref. [TU], is now performed. This is 
achieved here by later injecting a subsequent pulse with 
the local oscillator phase switched to measure the optical 
phase quadrature as in Fig. [lje) . Repeating this proto- 
col many times and post-selecting the measurement out- 
comes AQx within the desired window provides a pow- 
erful experimental platform to generate and fully recon- 
struct a non-Gaussian state of motion of a mechanical 
resonator. 

In order to prepare mechanical superposition states 
with T x there needs to be sufficient mechanical displace- 
ment induced optical rotation such that two distinct posi- 
tions give the same amplitude quadrature outcome. This 
is best achieved if the mechanical mean position gives 
zero rotation. For mechanical states that have a non- 
zero mean, which could have been conditionally prepared 
with a prior pulse |10j . non-Gaussian state preparation 
and tomography can be performed by providing a feed- 
back phase-shift (indicated by the arrow in Fig. life)) to 
rotate the optical scimitar to be centered about the Qx 
axis, as m Fig.fljd). Additionally, it is remarked that for 
optical rotation beyond that considered in ^ existing ex- 
perimental calibration procedures and the interpretation 
of optical phase measurements will require modification 
to take the optomechanical non-linearity into account. 

Studying the decoherence of quantum superposition in 
a mechanical resonator is an important step to deter- 
mine the feasibility of optomechanical systems as com- 
ponents for quantum information applications. Propos- 
als for such applications are numerous and include: quan- 
tum memory [24] , optomechanically mediated qubit-light 
transduction [25] and coherent optical wavelength con- 
version [26], to name a few. There is much literature on 
the topic of environmental coupling and decoherence [27] 
and so no detailed discussion will be provided here. How- 
ever, in the context of this proposal, what is important 
is the parameter n/Q, where Q is the mechanical quality 
factor. This parameter quantifies the rate of rethermal- 
ization normalized to the mechanical frequency and must 
be much less than unity to study the evolution of quan- 
tum mechanical phenomena over the timescale uf^ . For 
a temperature of 25 mK accessible with dilution refrig- 
eration and a Q — 5 x 10 6 this gives n/Q = 0.05 using 



the mechanical frequency above. With the full quantum 
state tomography available here this scheme allows the 
dynamics of mechanical superposition states to be mea- 
sured, which may be used to characterize the couplings 
responsible for decoherence thus allowing for improved 
mechanical device engineering. 

Furthermore, the significant mass involved in the spa- 
tial superposition offers a parameter regime that allows 
for an experimental test of collapse models. Very recent 
proposals in matter- wave interferometry [141 128j , which 
also consider the use of filtering type measurements to 
generate superposition, may provide the ability to test 
continuous spontaneous localization (CSL) [29]. The me- 
chanical resonator parameters considered here are not 
suitable to test CSL, predominantly as the superposition 
separation is small 30 . However, the separation can be 
larger than the distribution of the mass contained within 
the nucleus and so this can be used to test gravitational 
collapse [31 j - For example, using the parameters above 
(S = 2.0, xq = 10 fm) the separation is approximately 
28 fm and the diameter of a 28 Si nucleus is approximately 
8 fm. It may be useful in such an investigation to start 
with an initial squeezed mechanical state, as is consid- 
ered in Fig. [2jg-i) , as one can study a larger range of 
superposition separations as the probability density of 
measurement outcomes is broader. 



V. CONCLUSION 

This work has provided a means to measure displace- 
ment or displacement squared of a mechanical resonator 
using the optomechanical interaction linearly propor- 
tional to mechanical displacement by simply changing 
the phase in optical homodyne measurement. Displace- 
ment squared measurements have so far been predomi- 
nantly considered in dispersive optomechanics, however, 
the optimal square displacement measurement strength 
obtained in the scheme introduced here scales more fa- 
vorably than that available in dispersive optomechanics. 
This opens the possibility that optomechanics with a lin- 
ear mechanical coupling may also provide a route to ob- 
serve mechanical energy quantization, as was considered 
in [2]. Moreover, as was proposed in [15], with an X\ ;l 
coupling to a mechanical resonator one can prepare a 
superposition of positions via measurement. This, ap- 
plied to the X]^ coupling achieved here and combined 
with the ability to perform mechanical state tomography 
with time [TU] provides an alternative to Refs. [T7] and 
[19) to generate superposition of a mechanical resonator 
without the need for large single photon mechanical dis- 
placement gy m /K. Such mechanical superposition states 
are important to experimentally investigate in order to 
determine the feasibility of mechanical resonators as el- 
ements in quantum information applications and to ex- 
plore decoherence mechanisms arising from environment 
interaction or, for example, gravitationally induced col- 
lapse. 
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